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We investigate the possibility of using a dissipative process to prepare a quantum system in 
a desired state. We derive for any multipartite pure state a dissipative process for which this 
state is the unique stationary state and solve the corresponding master equation analytically. For 
certain states, like the Cluster states, we use this process to show that the jump operators can be 
chosen quasi-locally, i.e. they act non-trivially only on a few, neighboring qubits. Furthermore, 
the relaxation time of this dissipative process is independent of the number of subsystems. We 
demonstrate the general formalism by considering arbitrary MPS-PEPS states. In particular, we 
show that the ground state of the AKLT-model can be prepared employing a quasi-local dissipative 
process. 



I. INTRODUCTION 

Preparation of entangled pure quantum states is of in- 
terest in the context of both quantum information and 
condensed matter physics. In quantum information en- 
tangled states of qubits can act as a resource for quan- 
tum computing, e.g. Cluster states [l[ in measurement 
based quantum computing Q, while in condensed mat- 
ter physics entangled states represent ground states of 
strongly correlated systems. 

A possible scenario for the preparation of entangled 
states of interest is cooling the system to the ground 
state of an appropriate many body Hamiltonian. Alter- 
natively, we can generate a state of interest from an ini- 
tial pure state, e.g. a product state, which can easily be 
prepared with available resources. This is achieved either 
by coherent evolution generated by a system Hamiltonian 
(i.e. a sequence of quantum gates), or, more generally, by 
applying the most general physical transformation, which 
is mathematically represented by a completely positive 
map. We will discuss here another route: the preparation 
of entangled states by designing dissipative processes, so 
that we drive the system via non-equilibrium dynamics to 
a pure entangled state of interest \*f?) for long times, for 
any initial mixed state represented by a system density 
operator p, i.e. 

pi=^|*><*|. (i) 

In particular, we will consider a situation where the time 
evolution of the system coupled to a reservoir can be 
described as a Quantum Markov Process with dynamics 
obeying a master equation, 

P = C{p) (2) 
= -i[H, pl + ^ge hc e pc\ - c\cip - pc\c^J . 

£ 

Here, H represents a system Hamiltonian while the Li- 
ouvillian, L(p) = J2e 9i ( 2cipc\ - c\cgp - pc\cA , can al- 
ways be written in Lindblad form with ci a set of "quan- 
tum jump operators" and dissipation rates gi > 



Such a description in terms of a master equation is valid 
provided the system dynamics is slow on the time scale 
of the reservoir correlation time, as is the case for typical 
quantum optical systems. 

Thus the goal of dissipative entangled state preparation 
is to design quantum reservoirs and system-reservoir cou- 
plings, and to identify necessary and sufficient conditions 
for the master equation |(3J), such that the desired pure 
state of a many body system is obtained as the unique 
stationary state. Indeed, it will be shown below that 
for any given there is a master equation which yields 
the required state as the unique pure state within a relax- 
ation time T re i ax ~ 1/ min g^ , independent of the number 
of qubits. 

We will be particularly interested in a situation where 
qubits (or spin-1/2 particles) reside on a lattice. Thus it 
is natural to restrict dissipation represented by the Liou- 
villian L to quasi-local jump operators acting only on a 
neighborhood of a given qubit, which raises the question 
of the class of states which can be prepared with these 
resources. We will show below that examples of states 
which can be generated include stabilizer states, matrix- 
product states (MPS) or projected entangled pair states 
(PEPS). 

The results of the present paper are also of immediate 
relevance for a non-equilibrium condensed matter physics 
where pure many body states and quantum phases are 
prepared as a result of a driven dissipative system dy- 
namics. In a standard equilibrium situation of condensed 
matter and cold atom physics, states close to the ground 
state of the Hamiltonian, H \G) = Eq \G), are prepared 
by cooling the system, where in particular for a finite 
system p ~ e - H / k s T _ > |Q) (G| for temperature T — > 0. 
This has already lead to the preparation of intriguing 
quantum phases [J, [H, H, 0, H, H, In contrast, we 

obtain a pure state representing a non-equilibrium quan- 
tum phase as a result of the dynamics |T]) with the mas- 
ter equation (J2J). In recent work we have provided ex- 
amples of master equations, including the example of 
non-interacting bosons and paired interacting fermions 
corresponding to cold atoms moving in an optical lattice 
which are driven by coupling to quasi-local reservoirs into 
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pure states exhibiting long range order. While the goal 
of Ref. [ll[ was to study non-equilibrium condensates, 
Luttinger liquids and Kosterlitz-Thouless phases, in the 
present work we will provide the uniqueness proofs for the 
corresponding driven dissipative dynamics. Furthermore, 
the results of the present paper are of direct relevance for 
non-equilibrium spin models. As an example we will dis- 
cuss a master equation whose unique steady state is the 
ground state of the familiar AKLT model [12j |. 

The outline of the paper is as follows. First of all, 
we summarize some properties of the master equation. 
In Section III we characterize all stationary pure states. 
Since we are interested in unique stationary states, we 
derive a sufficient condition for the uniqueness of the 
stationary state. In Section IV we derive, for any mul- 
tipartite state, a dissipative process which can be used 
to prepare this state. That is, we construct a dissipative 
process for which the desired pure state is the unique 
stationary state. For this process it is simple to solve 
the underlying master equation analytically and to show 
that the relaxation time of the system is independent 
of the number of subsystems. In Section V we finally 
show that for certain states, like the 2D-Cluster states 
[l|, this construction can be used to choose jump op- 
erators quasi-locally. We furthermore derive a quasi- 
local dissipative process which has a general PEPS [l3[ 
as the unique stationary pure state. Also in the context 
of PEPS, we consider the ground state of the familiar 
AKLT-model [l2j and derive the dissipative process for 
which this state is the unique stationary state. We also 
prove the uniqueness of the driven noninteracting BEC 
and the 77-condensate of paired fermions given specific 
dissipative processes, complementing the work done on 
these states in fllj ]. 

II. PRELIMINARIES AND NOTATION 

We are interested in the stationary solutions of the 
master equation p = C(p), presented in Eq. @. Wc 
write C{p) as C(p) — £{p) — Q^p — pQ, where £ (p) = 
2^2 l gi2cipc\ is a completely positive map and Q — 

P — iH, with P = gicjci a positive semidefinite op- 
erator. Sometimes we denote £, as given in Eq. 
by C^H,ci} an d by £{ Ci } if we consider a purely dissipa- 
tive process. Note that the partition in the Hamiltonian 
and the dissipative part is unique if the operators q are 
traceless and orthonormal [Til ]. 

Since the master equation is linear, the eigenvalue 
equation C(a) = Act can be written as a matrix equation. 
Due to the fact that tr(£(cr)) = for any a, the eigenvec- 
tors to eigenvalues different than zero must be traceless. 
The eigenvalues can be complex, however, the real part 
of the eigenvalues is not positive (see for instance [la]). 
Considering a purely dissipative process, with hermitian 
quantum jump operators c;, all the eigenvalues are real. 
This is due to the fact that in this case the matrix cor- 
responding to the Liouvillian is hermitian. However, in 



general the matrix is not hermitian, it is not even diag- 
onalizable and has therefore generalized eigenvectors, a 
Furthermore, since C(a^) = (£(cr))t the eigenvalues 
occur in pairs of the form A, A* with the corresponding 
eigenvectors a and . The set of all proper and all gen- 
eralized eigenvectors, {<7i}, forms a basis in the operator 
space. Expanding p(0) in this basis, p(0) = ^ CiOi we 

obtain p(t) — J2iJ2j c l(^) cr j eXit where c^(i) are polyno- 
mials of degree less than the largest order of the Jordan 
block corresponding to eigenvalue A^ [1 €31 ] - 

We are interested in the stationary states of the evo- 
lution. That is we want to find the states, p, for which 
£{H,c,}(p) = 0. In order to do so, we use the following 
notation. By K(X) \R(X)] we denote the kernel [range] 
of a hermitian operator X and r{X) denotes the rank of 
X. Note that the pure states, which are not affected at 
all by the dissipative process coincide with the kernel of 
P, K{P) = {|*) such that cj |*> = 0V2 e L}. Thus, if 
I*) € D {H,ci} = K(P)C\ES(H), where ES(H) denotes 
the eigenspace of H, then C{ HjCl y(\^) (*|) = 0. This 
implies that any state p with R(p) = span{|0j)} such 
that C Df H ci y is a stationary state. We call these 

states dark states. Whenever it is clear from the context 
we omit the subscripts {H, c;} and write for instance sim- 
ply D for the subset of dark states. 

III. PURE STATIONARY STATES AND 
UNIQUENESS OF STATIONARY STATES 

In this section we first of all characterize all pure sta- 
tionary states. We show that a pure state is a stationary 
state of some dissipative process iff it is a dark state of 
some other. Then we derive a sufficient condition for the 
uniqueness of stationary states. 

Let us first of all, consider a given dissipative process 
and characterize the pure states, which are stationary 
states of it. 

Theorem 1. Let C be defined as in Eq. (0). Then 
£(1$) ($|) = iff the following two conditions are ful- 
filled: 

(1) \ X ) = Q ] |$) = A|$) for some A ed . 

(2) a |$) = A; |$) V; S L. for some A ; e <X with 
^2i9i\^i\ 2 = -Re(A), where Re(x) denotes the real 
part of x. 

Proof. ($|) = iff 

2^> |*,) = ffl*) ($|) = |$) ( X \ + \ X ) <<&| , (3) 

where \^i) = q 1$). Therefore, the operator A = 
\&) (xl + Ix) (^1 must be positive semidefinite. It can be 
easily verified that A > iff r{A) < 1. Thus, \\) = A |$) 
for some A e <T and A = 2Re(X) |$> ($|. The fact that all 
\^i) are in the range of £(\&) ($|) implies then, that Eq. 
© is fulfilled iff = A ; |$) with Y,i9i\M 2 = #e(A). 

□ 
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Using this theorem it is now easy to characterize all 
pure states for which there exists a Liouvillian such that 
the state is a stationary state of the master equation. 
Defining c\ = c\ - hi and H' = H — i^i 9ih{ c 'i) f + 
i^2i9i^i c 'i it i s straightforward to show that the condi- 
tions above imply that 1)H' |$) = A |$), with X e R and 
2)c\ |$) = VZ G {1, . . .m}. Thus, we have that there 
exists a Liouvillian C such that £(|$) ($1) — iff there 
exists a set of operators {cj}; = i.... m , a Hamiltonian H' 
and some A G R, such that the following two conditions 
are fulfilled: 

(1') c£|$) = 0VZ€l,...m 

(2') #'|$) = A|$). 

Note that the two conditions (1') and (2') are equivalent 
to |$) € DiH'.c'j- Therefore, |<&) is a stationary state, iff 
it is a dark state for some other physical process. Thus, 
in order to design a dissipative process, which leads to the 
desired pure state, we have to find a set of operators (or a 
single operator) which have only one common cigenstate, 
the corresponding eigenvalue can be chosen to be zero. 
Due to the results presented above, we know that the 
corresponding dissipative process will have the desired 
state as the unique pure stationary state. However, since 
we want to use this process for state preparation, we have 
to guarantee that there exists no mixed stationary state. 
How this can be ensured will be shown next. 

We consider the general master equation given in Eq. 
@ • Note that the corresponding set of dark states is in 
general a set, not a subspace. However, since we want 
to use the dissipative process for state preparation, or to 
drive the system to a certain (higher-dimensional) sub- 
space, we will consider here the situation where D is a 
subspace, i.e. all the states in D correspond to the same 
eigenvalue of the Hamiltonian. We are going to show that 
if there exists a stationary state, p, which is not a dark 
state, i.e. R(p) D, then there must exist a subspace 
of the Hilbert space, Ti, which is left invariant under the 
operators {q}. 

Theorem 2. If there exists no subspace S C Ti with 
S _L D such that ciS C S Vq, then the only stationary 
states are the dark states. 

Proof. We prove the statement by contradiction. That 
is we assume that there exists a state p, with R{p) = 
S' D such that C{p) — and show that this implies 
that there exists a subspace S C Ti with S _L D such 
that ciS C S Vq. Using the notation of Eq. (J2J) and the 
fact that £ is a completely positive map, we have that 
£(p) ~ Q^P + PQ ^ 0- Since £{p) is a positive semidef- 
inite operator we have that R{£{p)) = spanjq |£) ,1 £ 
L, |£) G S'}. We are going to show now that R(£(p)) 
must be within S' which implies that qS" C S' Vq. Since 
Q'p + pQ must be positive semidefinite K{Q} p + pQ) — 
{|$) such that ($|<?V + PQ\®) = 0} 2 K(p)(jD. 
Therefore we have that K(£(p)) D K(p) [JD which im- 
plies that R(£{p)) C R(p)f]D ± , where D 1 denotes the 



orthogonal complement of D. We have shown now that 
{q |0 ,Z e L, \£) £ S'} C S'f)D ± . It remains to show 
that if such a S' D exits, then there exists a subspace 
S C fulfilling that qS C S Vq. Note that the set 
D' = {|*) such that Qt |$) = j w i t h A e R} is 

equal to D [35|. Using again that Q^p+pQ must be posi- 
tive semidefinite which implies that ($| Q< p + pQ |$) = 
iff (Qt p + p Q) |$) = o we have that QV |$) = iXp |$) 
V |$) S D. Thus, p |$) S D V |$) S D, which implies that 
p can be decomposed as p = Pd + Pd ± i where R{pd) C D 
and R(p D ±) C D x . Now, since £(p) = £(p D ±) we only 
have to chose 5 = R(p D ±). 

□ 

This proof shows that if there exists a stationary state, 
p which is not a dark state, i.e. R(p) ^ D then p = 
Pd + Pd 1 - with R(pd) Q D and R(p D ±) C D 1 - such that 
R(Pd ± ) is invariant under the operators {q}. 

In the next sections we use the results presented above 
to design dissipative processes which can be used for state 
preparation. That is, we derive the jump operators, such 
that the system is driven into the unique stationary state. 
Due to the results above this state can be chosen to be a 
dark state of the process. 

Note that this goal can never be achieved using her- 
mitian jump operators q. The reason for that is the fol- 
lowing. If all operators q are hermitian, then the master 
equation can be written as 

£(p) = ^2-i[H >P ]+gi[[ci,p} ) c l }. (4) 

Thus, any operator which commutes with {q} and H is a 
stationary state, for instance, the completely mixed state 
cx 1 is stationary. 

IV. CONSTRUCTION OF DISSIPATIVE 
PROCESSES LEADING TO THE DESIRED 
STATE 

In this section we show that one can design the system- 
reservoir coupling in such a way that any multipartite 
state can be obtained as the unique stationary state of 
a dissipative process. Furthermore, we solve the corre- 
sponding master equation analytically by deriving the 
whole spectrum of the corresponding Liouvillian. Given 
this solution it is then easy to compute, any relevant 
quantity of the process, like for instance the relaxation 
time, or any correlation function. In fact, we show that 
the relaxation time is independent of the number of sub- 
systems. Even though these facts are not very surpris- 
ing, we will use the constructed process to show that for 
certain states the jump operators can be chosen quasi- 
locally. In general such a construction will not lead to 
a quasi-local dissipative process. Therefore, we demon- 
strate how the quasi-local operators can be constructed 
given a quasi-local description of the state at the end of 
this section. We start out by considering n-qubit states 
and generalize later the formalism to a d-level system. 
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For any n-qubit state 1^) = U\Q) n we construct 
a set of operators {cfe}JJ =1 such that the unique sta- 
tionary state of the dissipative process, described by 
p = £ {ck] (p), is |*). Note that |0>® n = |0 n ) is the unique 
stationary state of the dissipative process with jump op- 
erators dk = <J_ = li <g) (g) |0) fc (1| (g) l n , i-e. dk acts 

non-trivially only on the fc-th qubit. This can be easily 
seen, since dk \4>) = Vfc iff \<p) = |0 n ), which implies that 
|0") is the unique pure stationary state (see Theorem[T]). 
Furthermore, it is straightforward to show that for any 
|x) -L |0 n ) there exists a monomial of the jump opera- 
tors, P(|x)) such that (0"| P(|x)) |x) ^ 0. This shows 
that |0 ra ) is the unique stationary state of C{d k }- 

We construct now the operators Ck which lead to the 
unique stationary state \*f?) = U\0 n ). Defining Ck = 
Ud k W we have c k \<j>) = iff \<j>) = U |0"> = |#). One can 
use the same arguments as above to show that this is the 
unique stationary state. This immediately implies that 
if the jump operators Ck can be written as U(li t ,,,k-i <g> 
|0) fc (l| <g> l k+ i tn )W and I*) = U\0 n ) is an eigenstate 
of the Hamiltonian H, then is the unique stationary 
state of L{H,c k }- 

Note that the statements above would not be changed 
if we would use an invertible matrix, X, instead of a 
unitary, as long as this does not lead to more common 
eigenstates of the operators {XckX^ 1 }. 

A. Analytic Solution of the Master Equation 

In this subsection we solve analytically the master 
equation p = £{ Ck y(p), for any set of operators {ck} 

which are unitarily equivalent to the set {o_ }, i.e. 

Ck = Ucr_U^ for some unitary U. 

In order to compute the eigenvalues and eigenvectors 
of the Liouvillian £{ Cfc } we note that C^ Ua _( k ) uf y(a) — 

U£r W JWaU)W. Therefore, £ {£/o . W[/+} (ct) = Act iff 
C^h)^{WaU) = \W<jU. Thus, computing the eigen- 
vectors and eigenvalues of Cr(k)-. gives us immediately 

the eigenvalues and eigenvectors of £{ Ck }- 

Since C^ k )^ describes the situation where n qubits 

are interacting with identical and independent reservoirs, 
we only have to find the eigenvectors and eigenvalues 
of C r (fc), for one fixed k. This describes the situation 

where a single two-level system is interacting with a 
thermal bath. The solution to the eigenvalue problem 
is well-known [13, (3 ■ We find £ {o . (fc)} (of } ) = A* of iff 

of £ E* = {of = |0) (0| ,of = of }|Jspan{of = 

of jof = Oy } and the corresponding eigenvalues 
are Ao = 0,A3 = —2gk, and the two-fold degener- 
ate eigenvalue is Ai = A2 = —gu- Thus, a basis 
of the eigenvectors of the total Liouvillian is {of ® 

■ ■■o^ with of e ^k}i k e{o,l,2,3} and tne correspond- 
ing eigenvalues are At t ,...,»„ = X^fe^'fc- ^°t e that the 



eigenvectors of £{ Ck } are just the unitary transformation 
of the eigenvectors of £,(*,), with the same eigenval- 
ues. Thus, writing the initial state in the eigenbasis, 

P(0) = UY, iu ... in C h,-in4l ] ® • Ul We 0btam 

p{t) = e Cf (p(0)) (5) 
= U X: c iu .., n e-^ ^>...CTf[/t. 

Let us stress here the fact that there exists no purely 
imaginary eigenvalue. If such a pair (recall that complex 
eigenvalues occur in pairs, see introduction) existed then 
the additional condition that lim^oo p(t) = p ss , where 
p ss denotes the stationary state, would not be satisfied. 
In order to see how fast the system is driven into the 
stationary state, we compute the relaxation time, T re i ax 
which is defined as the inverse of the minimum of the 
negative real part of A^....^ different than zero. We find 
Treiax = 1/rnirXfeOfe and is therefore given by the minimal 
coupling constant. Note that it is not very surprising that 
Treiax > l/mirifefffc) since the system cannot be driven 
faster into the stationary state. The reason why this 
amount of time is already sufficient is because the evolu- 
tion of the qubits can be decoupled. In other words, for 
fixed coupling constants, the relaxation time is constant 
in the number of qubits. 

It is straightforward to generalize this formalism to d— 
level systems. One simply has to replace the operators 
|0> (1| by the a matrix J d = |0) (1| + |1) (2| + |2) (3| + 
. . . \ d — 1) (<Z|, i.e. a Jordan matrix with only one eigen- 
state (here to eigenvalue zero) [36]. Only the state |0) is a 
proper eigenstate of Jd- All the other computational ba- 
sis states are generalized eigenvectors, which means that 
for any k ^ there exists an i such that J l d \k) = |0). 
This implies that for any state \<j>) there exists a poly- 
nomial of the operators Jd, P(Jd), such that the overlap 
(0| P{Jd) \ <f>) is not vanishing, which is exactly the prop- 
erty that we need in order to prove that the stationary 
state |0 n ) or respectively U |0 n ) is the unique stationary 
state. Since also in this case, the master equation decou- 
ples for the different subsystems it is straightforward to 
solve it analytically. 



B. Quasi— local dissipative Processes 

If one wants to prepare the state |^) using a dissipative 
process, then one will be interested in a simple physical 
interaction between the reservoir and the system. One 
requirement, for instance, could be that the operators Ck 
are quasi-local, which means that they act non-trivially 
only on a small number of qubits. Depending on the state 
|\&) the operators Ck might be chosen quasi-local, as we 
are going to show next. Note however, that this cannot 
be true for any state. The reason for this is the following. 

As shown before, for any n-qubit state |^) one can find 
n operators Ck which uniquely define the state |$), in the 
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sense that \*S?) is the only state which is a (right) eigen- 
state to eigenvalue of all operators c^. Thus, a descrip- 
tion of the state \^f) is the set of operators {efe}£ =1 . In 
general one might need more than one operator per sub- 
system. We denote the corresponding set of operators 
by {c£}, where a = l,...,d, for some d, and call this 
set a quasi-local description of if all c£ are quasi- 
local. The Quantum Kolmogorov complexity, i.e. the 
number of classical bits required to describe the state 
|\&) equals the classical Kolmogorov Complexity of the 
set {c£}. If all these operators would be quasi-local, 
then the Kolmogorov Complexity scales only polynomi- 
ally with the number of qubits. It is known however, 
that for any n there exists a n-qubit state, whose Quan- 
tum Kolmogorov Complexity scales exponentially with 
the number of qubits [19(. 

Before we discuss more generalized schemes we want 
to use the process discussed above to show that cer- 
tain states can be prepared using quasi-local dissipa- 
tive processes. We consider the state |W) = U\0) ™. 
Now, suppose that U = U\U?...U n , where each of the 
unitaries Uk commute with each other and are quasi- 
local. Let us assume that Uk is acting on particles 
k — l,k,k+l. Then the jump operators = U<?kU* = 
Uk-iUkUk+i&k{Uk-iUkUk+iy are also quasi-local and 
the dissipative process corresponding to Ck has as a 
unique stationary state the state \^f). As shown before, 
for such a process the relaxation time is constant, i.e. 
independent of the number of subsystems. 

For certain cases such a simple construction will not be 
possible. Therefore, we describe here a general method 
of deriving the jump operators, i.e. the dissipative pro- 
cess, which give rise to the desired state. Since the jump 
operators are not hermitian (at least not all of them can 
be chosen to be hermitian), they might not be diagonal- 
izablc. However, one can use the Jordan normal form 
to gain some insight in the necessary properties of these 
operators [20(. The jump operators must be chosen such 
that there exists only one common eigenstate to eigen- 
value zero. The Jordan decomposition of a matrix c is 

c^SJS- 1 , (6) 

where J denotes the Jordan matrix of c and S is an in- 
vertible matrix. J is a block diagonal matrix with di x di 
diagonal blocks 

/Xi 1 \ 
A l 1 

JM= ■■. ■■. ■ (?) 

Aj 1 

V \i J 

The number of Jordan blocks is the number of lin- 
early independent eigenvectors of c. Let us for sim- 
plicity consider here the case where S is a unitary, U. 
The proper eigenvectors of c are then f e/(k)/j where 

/(&) — Ei=o ^ + ^ ano - l efe ) denotes the standard ba- 
sis. A simple example is the operator er_. The only 



eigenstate is the state |0). Considering multipartite en- 
tangled states, the operators must have more than a 
single eigenstate to eigenvalue (like 1 ® ® 1 has). 
Therefore the matrices cj, will have Jordan matrices with 
several Jordan block and all eigenvalues 0. The unique 
eigenstate which is common to all operators is the state 
one wants to prepare, |\&). Apart from that, one has to 
assure that any other state can be mapped into some 
state having non- vanishing overlap with 

One might also chose a single jump operator, which has 
only one proper eigenstate (which corresponds to eigen- 
value zero) (see also (l5j). In order to be more precise, 
let us denote by {l^i)}^^ 1 an orthonormal basis with 
\4>o) = \^) being the state we want to prepare. The op- 
erator C — X)i=o 1 \4>ii) {4>i+i I has only one proper eigen- 
state, namely \^>). For the generalized eigenvectors, \<j>i), 
for i > it holds that C { |<fo) = |*). Thus, for any 
state |x), there exists a k such that C k \x) y^ 1 0. This 
shows that the single operator C leads, as the operators 
{c^} described above, to the unique stationary state 
However, in contrast to the operators {ck} the operator 
C will act non-trivially on all subsystems and is therefore 
(in general) harder to implement experimentally. 

In summary, a sufficient condition for the existence 
of a quasi-local dissipative process, which can be em- 
ployed to prepare a certain state is the following. First 
of all, the state must have a quasi-local description, i.e. 
there exists a set of operators, {c£} such that the only 
common eigenstate is \^f). This implies that there ex- 
ists only one dark state for the corresponding dissipa- 
tive process. Second, if there exist polynomials Pi{{c^}) 
such that the states -P/({c£}) form a basis in the 
Hilbert space, then \^f) is the unique stationary state. 
More generally, if there exist polynomials Pi({c%}) such 

that smm = e^(K})i*)<*i^(K}).= x 

where X is invertible, then \*Sf) is the unique stationary 
state. This can be seen as follows. If £ (|*) (*|) = X, 
with X invertible, then we have for any state \x) € H, 
(X| £(\#) (I'D |x) # 0. Since (*|) is a sum of posi- 

tive semi-definite operators this implies that there exists 

at least one term, P/({c£}) l*> (*l p i({ c k})> which has 
a non-vanishing overlap with the state |%) (x|- 



V. EXAMPLES 

In this section we illustrate the formalism described 
above by applying it to several examples. In subsec- 
tions [V2] and [VD] we show that the dissipative processes 
which we analyzed in detail in [Tlj . have unique station- 
ary states, namely the BEC-state and the ^-condensate 
respectively. In subsection IV Bl we derive the processes 
for stabilizer states, and in subsection IV CI we consider 
arbitrary PEPS states and show how, for instance the 
ground state of the AKLT-model can be generated with 
a quasi-local dissipative process. 
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A. Example: Driven Dissipative Hubbard 
Dynamics of Bosons on a Lattice 

In Ref. (lH | we have described a driven dissipative Hub- 
bard dynamics of bosonic particles on a lattice. The cor- 
responding dynamics was written in terms of a master 
equation ([2|) with H a Hubbard Hamiltonian, contain- 
ing the coherent hopping of particles between the sites 
of the lattice, and their interaction, and where a Liou- 
villian with quasi-local jump operators was designed to 
drive the system into a non-equilibrium condensate of 
bosons or paired fermions. While Ref. [ll[ focused on 
non-equilibrium condensed matter aspects and in partic- 
ular on the effect of interactions, and implementation of 
this master equation with cold atoms, we claimed - but 
did not prove - the uniqueness of the many body dark 
states for this master equation. This missing proof will 
be provided here based on the theorems derived in the 
previous sections. 

We consider N bosons hopping on a d-dimensional lat- 
tice with Hamiltonian 

H B = H + V = - J £ a\ aj + ±C/ ]T , (8) 

which is written as the sum of a kinetic and interaction 
energy term. Here Oj denotes the destruction operator for 
a boson on site i, and (i, j) denotes adjacent lattices sites. 
Ref. suggested that for non-interacting bosons (U = 
0) a master equation (f5]) with the above Hamiltonian and 
quasi-local dissipation 

ct = Cij = (a\ + af\ (a, - Oj) (9) 

has the unique steady state solution in form of a pure 
Bose Einstein condensate representing a state with long 
range order 

|BEC)=a q f /^|vac). 

Here a q = l/V M d Oj-e" 16 ^ is the destruction oper- 
ator for quasimomentum q on a lattice with M d lattice 
sites, lattice vectors e\ and spacing a. It is easy to see 
that the state |BEC) satisfies 

(i) V(i,j) K - Oi) |BEG) = 

(ii) H |BEC) = JVe q =o |BEC) 

with e q= o the single particle Bloch energy for quasimo- 
mentum q = 0, and thus is consistent with the conditions 
of a pure steady state of the master equation. We show 
now that this state is unique. 

Before proceeding we find it convenient to go to a 
quasimomentum presentation. The dissipative part of 
the master equation then takes on the form @ with jump 
operators 

ct = Cq , A = £(1 + e^-^Xl - e- ike >)4_ q a k 

(10) 



and couplings constants, gi independent of I. 

To show uniqueness we show first off all that the 
BEC-state is the unique dark state and that there ex- 
ists no invariant subspace for our choice of jump oper- 
ators. For a fixed particle number N, the first term 
in a is a creation operator and has no eigenvalues; in 
particular no zero eigenvalues. Thus, in order to iden- 
tify dark states \D) with zero eigenvalue, we may re- 
strict ourselves to the equation (a* — a,j)\D) = V(i,j). 
Taking the Fourier transform, this translates to (1 — 
e"^)a q |L>) = Vq. Thus the BEC state with q = 
is the only dark state. Next, we will construct for ev- 
ery state |$) in the Hilbert space a polynomial opera- 
tor 0(cq t \), where the jump operators c q; A are given in 
Eq. (JTDJll such that (BEC| 0|$) ^ 0. With the nota- 
tion n = {. . . , rik, . . .}, the states |n) = n q ( a q)" q l u ) 
form a basis in the Hilbert space, a general state can be 
written as |$) = J2fn\ n )- We select a state |m) with 
/m 7^ and the number of particles in the zero momen- 
tum mode maximal. Furthermore, for each q =/= we fix 
A q such that eA q is not orthogonal on q, i.e., eA q q ^ 0. 

Note that (BEC|(c q , Aq )"i(aj c ) mk |0) ^ only if n q = m k 
and q = k. Thus, applying the polynomial operator 

= n q ^ (c q ,A q ) q to |$) provides a finite overlap with 
the BEC, which implies that the BEC is the unique sta- 
tionary state of the dissipative process. 

B. Preparation of Stabilizer States 

As a second example we show that stabilizer states can 
be obtained as the unique stationary states of a dissipa- 
tive process involving only quasi-local interactions. 

We denote by X, Y, Z the standard Pauli opera- 
tors. The Pauli group, V, consists of all Pauli matri- 
ces, X, y, Z, 1, together with the multiplicative factors, 
±l,±i. V n = "P®" defines the Pauli group on n qubits. 
A n-qubit state, \^>), is called stabilizer state if it is 
uniquely defined as the only eigenstate to eigenvalue +1 
of a hermitian subgroup (of order n) of Vn, called the 
stabilizer of 1*) [Hi]. That is, |^) is the only state left 
invariant under the subgroup 

S(|*» = {5GP„: 5 |*) = |*». (11) 

Let us denote the generators of this group by Ui, 

1 = 1, ...n. A subset of the stabilizer states are the 
so-called Graph states [22J, which are associated to a 
mathematical graph which consists out of vertices and 
edges. Whenever two vertices a, b are connected by an 
edge, we say that b is in the neighborhood of a and write 
b G N a . A n-qubits Graph state, |^ r o,...,o) 3 can now be 
defined as the unique eigenstate of a set of independent 
commuting observables C4 = X k Z Nk = X k II&eAr fc z b, 
where Wf, — , for W € X, Y, Z denotes the Pauli 
operator a w acting on qubit k. Note that these uni- 
taries define a unique basis, the so-called Graph state 
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basis, which we denote by {|^ , i 1 ,...,i„)} Jj e{o,i}- It can be 
shown that |* 41 ,..., 2 „) = Zj 1 ® . . . Z% |* ,...,0> i 22 ], where 
Uk\^i lt ...,i n ) = (— l) 1 * l^n, ...,»„)• In particular we have 
that 14^,. ..^=0,...^) = Z l k k 1^,...,^=!,...,^). We con- 
sider now those Graph states for which all the unitaries 
Uf. are quasi-local. An example would be the linear Clus- 
ter state, where Uk = Zk-\XkZk+\- An other example 
would be the 2D-Cluster state, which is a universal re- 
source for quantum computations Q. As shown in Q, 
once such a state is prepared any quantum computation 
can be performed by means of local measurements only. 

We define the operators Cfc = ^(1 + Uk)Zk, which 
act only non-trivially on the same qubits as Uk and 
are therefore quasi-local. Note that 1 + Uk is a projec- 
tor onto the subspace of eigenvectors of Uk with eigen- 
value +1 and Zk changes any eigenstate | v E'i lj ...,i fc =o,...,i„) 
to |^i 1 ,...,i fc =i,...,i n ) and visa versa. Thus, Ck is an op- 
erator which has only one 2™~ 1 -fold degenerate eigen- 
value, namely and the corresponding eigenspace is 
span{|*i 1 ,...,i fe=0 ,...,i„> , with ij e {0,1}, for j ± k} 
Thus, we have that Ck \ <j>) — iff Uk \(j>) = \4>), which 
implies that the only state which is an eigenstate to all 
operators Ck is the Graph state. Using the definition of Ck 
it is easy to show that Ck — Ua_ W , where U transforms 
the computational basis into the Graph basis. Due to the 
results presented in Sec. Ill, this shows that the Graph 
state is the unique stationary state of the process de- 
scribed by £{ Ck }- Since any stabilizer state is up to some 
local unitary (actually local Clifford) operation, V, equiv- 
alent to a Graph state, the operators c(S)k = VckV^ are 
quasi-local and C{ c (s) k } nas as unique stationary state 
the stabilizer state. Furthermore, the relaxation of this 
process in independent of n and there exists no purely 
imaginary eigenvalue of the Liouvillian. 

Note that we can write Ck — ZkPk, with Pk = ^(l— Uk) 
is the projector onto the eigenspace to eigenvalue —1 of 
Uk- Thus, the evolution corresponding to these operators 
can be implemented using a feedback mechanism. 

Another way to show that the stabilizer states can be 
prepared using quasi-local jump operators would be to 
use the relation between the unitary which generates the 
state and the operators {ck}, as described in Section IV. 



C. Preparation of MPS-PEPS states 

As mentioned above, one has to choose non-hermitian 
matrices to guarantee that the desired state is the unique 
stationary state. Since otherwise, at least the completely 
mixed state is a stationary state too. The aim of this 
section is to demonstrate the general construction of the 
jump operators. Therefore, we consider the so-called 
product entangled pair state (PEPS)[23|. In one dimen- 
sion these states are called matrix product states (MPS). 
Let us consider n d-dimensional systems. Any state de- 



scribing these systems can be written as 

d 

|*)= £ tr^UiJ-.-^li!,...^), (12) 

i\ ,...i n — l 

\k] 

where A\ are D x D matrix, with D being the bond 
dimension [24-]. Using this way of presenting the states 
and especially the generalization to 2D, has been proven 
to be very powerful to determine, for instance, the ground 
states of some Hamiltonians |23j . 

Similar to the procedure above we design now the in- 
teraction between the reservoir and the system, such that 
the desired MPS-PEPS is the unique stationary state. In 
[2f| it has been shown that for many MPS-PEPS one can 
construct a frustration free Hamiltonian which has this 
state as the unique ground state. These Hamiltonians 
are of the form 

# = 5>fc, (13) 

k 

where the hermitian operators hk are quasi-local. The 
ground state, |\E r ), of these Hamiltonians is uniquely de- 
fined by the equations 

fcfc|*) = A mill |*) Vfc, (14) 

where A m ; n denotes the minimal eigenvalue of hk- That 
is, the ground state of the Hamiltonian corresponds to 
the ground state of the quasi-local Hamiltonians. 

In order to obtain these states now as stationary states 
of a dissipative process we construct the operators Ck 
which have a unique common eigenstate, which is the 
desired MPS-PEPS. In order to do so we consider the 
ground states of the operators hk, D(hk) (note that this 
must be more than 1 -dimensional) and again construct 
the non-hermitian operators Ck whose eigenvectors span 
D(hk). Then, the only common eigenstate of the oper- 
ators Cfc is |\&). To ensure that there exists no mixed 
stationary state, one might need to consider more than 
one operator per site. 

Let us illustrate the general idea by considering as an 
example the ground state of the AKLT-model [12| • His- 
torically, this state occurred first in the context of con- 
densed matter physics. There it was shown to be the 
unique ground state of a Hciscnbcrg-likc Hamiltonian 
(see below) [l2j. Recently this state has attracted in- 
terest in quantum information theory, due to its useful 
properties for quantum communication [26| . The Hamil- 
tonian has the following form: 

H = y^fofc, with hk = SkSk+i + g (§ k §k+i) ,(15) 

k 

where the operators S a with a € {x, y, z} denote the 
spin-1 operators. The quasi-local Hamiltonians, hk act 
non-trivially on system k and k+1. As mentioned before, 
the ground state of the Hamiltonian, which we denote by 
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|\&), is the ground state of all quasi-local Hamiltonians, 
h k . 

Since the Hamiltonian is translational invariant we 
omit the index k as well as the identity operator when- 
ever it does not lead to any confusion. The operators 
hk have the eigenvalues 4/3, which is 5-fold degenerate 
and —2/3, which is four-fold degenerate. We denote by 
\^k) (l^fe)) an orthonormal basis of the subspace corre- 
sponding to the eigenvalue 4/3 (—2/3) respectively. The 
operators Ck (acting on system k and fc+1) can be defined 
as: 

Cfe = |*l)<$l| + |*2)<*2| + |*3><*3| (16) 
+ |¥ 4 ) <«4| + |$4> <$ B | • 

It is easy to see that the only states which are eigen- 
states of Cfc are states in the subspace spanned by 
({\^ , i- e - the ground state subspace of hk- Thus, 

the set of operators {ck} is sufficient to guarantee that 
there is no other pure stationary state. 

In order to show that there is no mixed state, we note 
that the operators Ck can be written as Ck = UkPk, where 
Uk is a unitary matrix and Pk = |(2/31 + hk) is the 
projector onto the eigenstates of hk with eigenvalue 4/3. 
The choice of the unitary is by no means unique and we 
will use this freedom to show that also in this case there is 
no mixed stationary state. Let us write c£ = U^Pk, then 
cl |$) = Vfc, a iff |$) = |*). For any state |x) different 
than |*), there exists a k such that Pk\x) 7^ 0- Without 
loss of generality we assume k = n. We construct now a 
completely positive map, £k{p) = S Q ( c fc)V c fc su ch that 

£(|*> = £n O £ n -l 0...0 ^(1*) (tf |) CX 1 2 ,3,...,„ ® 

P n and therefore, (x|£(|*) (*|) |x) + 0. Since all the 
terms occurring in £ (|*) (*|) are positive semidcfinite, 
this implies that there exists a monomial P — c"" . . . c" 1 
such that ( X | Ft |$) (^| P\ x ) = \ (*| p \ x ) | 2 ^ 0. Thus, 
the state is unique. We choose sufficiently many U a such 
that J2 a (UkYpUk = \% ® tv k ,k+i(p) for any state p 
|37j |. Using that tr k(Pk) oc 1, we find 

£ 2 (£i(|*) (*|) oc £ 2 (Pi ® tn, 2 (|*) (*|) oc (17) 

ll®P2(8)tri,2,3(|*) (*D- 

Continuing in this way we end up with l2.3,...,n ® Pni 
which shows that the ground state of the AKLT-modcl 
is unique. 

D. The driven 77-condensate 

The second example for a dissipatively driven state 
given m is the 77-condensate of paired fermions. First 
described by Yang [28j, the 77-condensate is an exact ex- 
cited eigenstate of the Fermi-Hubbard (FH) model for 
fermions with two internal states, "f, J., on a bipartite lat- 
tice in d dimensions with M d sites. Its properties derive 
from the fact that the FH-Hamiltonian 

H FH = -J Y, fy^v+uYflifuflifu (18) 



and the 77-creation operator 

^-^E^)/ T VL (19) 

obey the commutation relation \ti\Hfh\ = Urj'f, where 
S(l) in P^]) alternates between ±1 in a checkerboard pat- 
tern. The f^ ; operators create a fermion with spin a =|, 
I on site I and obey canonical anticommutation relations. 
In (|18[) . J denotes the tunneling rate and U the on-site 
inter-species interaction. 

From this observation, it follows that the state |rj, N) = 

(rf) |0) i s an exact excited eigenstate of Hpu, with en- 
ergy NU, irrespective of the boundary conditions. In 
the following, we assume periodic boundary conditions. 
In position space the 77-condensate can be understood 
as a superposition state of all possible vectors \D^- dk ^), 
which denote position basis states in Hn that have 
only double occupation, at positions dk, k G [1, . . . , N]: 
\V,N) = E {dk} S({d k })\D^), where S({d k }) = ±1 if 
the number of even dk is itself even (odd) on the even- 
indexed sites. A pair of fermions on the same site is called 
a doublon in the following. 

The 77-condensate, while never a ground state of the 
Hubbard-model, is interesting from the perspective of 
many-body physics, as it exhibits perfect superfluidity 
in any spatial dimension, without any approximations. 
77-pairing {{rprj) ^ 0) in the ground state has also been 
considered for doped negative U Hubbard models (29j 
and extended Hubbard models (see e.g. [30|). 

We define the jump operators 

4 1} = (vj-vl +1 )(m + m+i) (20) 

c, (2) = n^ifljfu+i +n u+1 fl l+1 f u 

4 3) = (fit + fh+i)(fu + W)(l - n lt i)(l - n Ll+1 ) 

c, (4) = {flj 1,1+1 + fli +1 fl,i)nun u+ i 

on the Hilbert space Hat, in which all states have N spin- 
up and N spin-down fermions. Here, rjj = ft if} ;j an d 

^a. J, — f a ifa,l- 

A dark state |*) G TLn is defined through the condition 
c^ fc ^ I = V k, I. In the orthonormal position basis, I*) 
is defined through coefficients ^i 1 ,j 1 ;...-,i M ,jm i where ii = 1 
(= 0) if there is (no) spin- up fermion on I, and ji be- 
ing analogously defined for spin-down fermions. Besides 

normalization, all ^4's must obey 53/=i */ = Y^fLi 3l = 
N. In this notation, cj imposes that the coefficients 
for any dark state must obey -Ai.i ; o,o = — A),0;i,i an d 

^4-l,l;l,0 = ^1,1;0,1 = j4l,0;l,l = j4o,1;1,1 = V 7fc, jk, 

k 7^ 1,1 + 1 (the ife's and jk's in A are suppressed for 
brevity). The first of the two conditions is the essen- 
tial property of the 77-condensate: \D^ dk ^) differing in 

the position of one doublon by one site carry opposite 
(2) 

sign. c\ imposes the constraint Ai.o : o,i = -Ao.i:i,o = 0, 
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which signifies that fermions of opposite spin on adja- 
cent sites are associated into doublons (on the site of the 

spin-up fermion). c\ then imposes Ai i0; o,o = — A),0;i,o 
and j4i,o;1,o = 0. These jump operators induce a diffu- 
sion process, which delocalises the spin-up fermions pro- 
vided they do not encounter spin-down fermions or dou- 
blons. Furthermore two spin-up fermions may not sit on 
adjacent sites in these dark states, cf^ finally imposes 
^4i,i;i.o = — -Ai,0;i,ij i- e - spin-down fermions must also 
be delocalised over the lattice, but only atop the spin-up 
fermions. 

Using the conditions imposed by c^- 1 ', and c^ 3 -* 
and the constraint on the number of spin-up and spin- 
down fermions, it is straightforward to show that the 
^-condensate is the only dark state. Note that \rj,N) is 
trivially also a dark state to . We assume that there 
are nonzero coefficients in |\&) for configurations with un- 
paired spin-up and spin-down fermions. The conditions 
for the dark state then immediately yield a contradiction. 
This is the case, as the condition A^ojo.o = — A),0;i,o 
means that any unpaired spin-up fermion may be shifted 
in position if no spin-up fermion, spin-down fermion or 
doublon is on the adjacent site. By assumption and the 
condition, all the coefficients to these configurations must 
also be nonzero. Shifting sufficiently often, one of these 
obstacles will eventually be encountered - at which point 
the other conditions will yield a contradiction. 

Thus, any dark state may only have nonzero coeffi- 
cients for the doublon basis states \D^ dk }). The first 
condition from the then immediately yields that all 
these coefficients must be equal to S({dk})- This is the 
case, because starting from any doublon configuration 
with nonzero amplitude, repeated application of this con- 
dition yields that any other possible doublon configura- 
tion on the lattice must also be nonzero, with its sign 
also obeying the rules for the 77-condensate. 

The proof of uniqueness again requires that for any 
l^) € Hn a monomial P({c| fe ^}), k = 1,...,4, can 
be constructed s.t. (rj,N\P\^f) 7^ 0. The proof pro- 
ceeds differently for two different cases: a) For ^ 
span{|Z?{ dfe })}, where \D^ dk ^) denotes a basis state car- 
rying only doublons on sites dk, k G [1, . . . , N], b) For 
I*) S span{|D dfc )}. 

a) The proof for the uniqueness of the dark state in- 
dicates that we can proceed analogous to the case of the 
BEC, but now in position space instead of momentum 
space: Pick one position basis state l^ -*) with unpaired 
fermions occurring in |\&) with nonzero amplitude. It is 
clear from the above proof that we can always find a se- 
quence of applications of c[ 3 ''s and c! [ s s.t. a spin-up 
and spin-down fermion are on adjacent site (we observe 

that c[ 4 ^ allows a doublon to swap positions with a spin- 
up fermion on an adjacent site). Application of c^ 1 then 
associates these into a doublon. It is straightforward to 
see that this sequence can be constructed s.t. it maps 
\4>^) into a particular \D^ dk ^), and any state orthogonal 
to to zero. 



b) We can proceed directly analogous to the proof of 
uniqueness for the ground state of the AKLT-model, by 
generalizing the jump operators to Cj ° = U"P\, 
with Pi the projector on the symmetric state of one dou- 
blon on sites I and I + 1. The unitaries, U" now must 
conserve particle number. We define completely positive 
maps £ 1 again as in IV CI Applying them in sequence 
on the ^-condensate yields a density matrix with nonva- 
nishing overlap to |\&), which proves the existence of a 
monomial P leading to finite overlap, as for the ground 
state of the AKLT-model. 



VI. CONCLUSION AND OUTLOOK 

In this paper we have shown how to use dissipative 
processes for state preparation. That dissipation, in con- 
junction with measurements , can be used to prepare pure 
states is known (see e.g. jU HH), where usually single- 
or few-particle states arc considered. Here, we demon- 
strated how dissipation can be employed to generate mul- 
tipartite states. For an arbitrary n c£-level state we con- 
structed a dissipative process (with constant relaxation 
time), consisting of n jump operators, which has this 
state as the unique stationary state. For certain states, 
like the Cluster states, we showed that this process can 
be implemented using only quasi-local operations. Apart 
from that we demonstrated how a quasi-local dissipa- 
tive process, which is suitable for state preparation, can 
be constructed for a given state, which has a quasi- 
local description. We illustrated this method considering 
the ground state of the AKLT-model. Furthermore, we 
showed that the quasi-local dissipative processes, which 
we considered in [ll|, have as unique stationary states 
the BEC-state, and the 77-condensate, respectively. 

The processes discussed here might be also used to 
gain some insight in the computational complexity of cer- 
tain problems, like for instance the so-called satisfiability 
problem [331 ]. There, the aim is to find out whether there 
exists a common solution to a set of Boolean equations. 
To be more precise, the problem is to determine a 71- 
bit string for which m = poly(n) given clauses, involving 
only 3 (2) variables for the 3-SAT (2-SAT) problem re- 
spectively, hold true. It is straightforward to construct 
a dissipative process such that any computational state, 
corresponding to the classical solution of the problem, is 
a stationary state of the process. Using the formalism 
we developed here, we are going to investigate the differ- 
ence between the NP-complete 3-SAT problem and the 
2-SAT problem, which can be solved (classically) in poly- 
nomial time. It might well be that looking at this prob- 
lem from this completely different point of view leads to 
a new insight. Furthermore, it would be interesting to 
establish a relationship between the relaxation time and 
the non locality of the jump operators, which cause the 
system to evolve to the desired state. Apart form that, 
the investigation of dissipative processes, which lead to 
higher dimensional dark state subspaces might lead to in- 
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teresting applications for quantum computation and stor- 
age of quantum information. 

In [34[ the authors investigated independently similar 
aspects of dissipative processes. It has been shown there, 
that these processes can also be employed for universal 
efficient quantum computation. 
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